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ABSTRACT 


The  deformation  response  of  polycrystalllne 
AA2024-T4  to  cyclic  loading  has  been  determined  by 
conventional  and  double-crystal  dlffractometry. 

These  studies  show  a  preferential  work¬ 
hardening  near  the  surface  compared  to  the  bulk. 

The  defect  distribution  as  a  function  of  depth  was 
also  investigated  as  a  function  of  the  fraction  of 
fatigue  life.  From  these  studies,  it  was  determined 
that  the  remaining  fatigue  life  could  be  predicted 
by  determining  the  average  ratio  6/8*  where  8  is 
the  linewidth  at  any  number  of  fatigue  cycles  and 
8*  is  the  critical  linewidth  at  fracture.  Current 
studies  have  shown  these  predictions  to  be  valid 
for  tests  where  the  stress  amplitude  is  held  con¬ 
stant  for  the  entire  test.  A  few  tests  have  also 
been  performed  where  the  loading  sequence  is  more 
complicated.  In  these  tests,  the  X-ray  predictive 
method  was  found  to  be  far  more  accurate  than  con¬ 
ventional  methods. 


ADMINISTRATIVE  INFORMATION 

This  investigation  is  part  of  an  in-house  research  program  at  the 
David  W.  Taylor  Naval  Ship  Research  and  Development  Center  (DTNSRDC) .  It 
was  conducted  under  Program  Element  61152N,  Task  Area  ZR  022-0101,  Work 
Unit  2082-004. 


INTRODUCTION 

Early  investigations  of  fatigue  failure  in  metals  provided  substan¬ 
tial  evidence  of  the  surface  sensitivity  of  mechanical  behavior.  Surface 
studies  revealed  the  formation  of  slip  bands  and  their  topological  devel¬ 
opment  into  intrusions  and  extrusions  associated  with  subsequent  crack 

1-4* 

initiation  and  eventual  failure.  Other  fatigue  investigations 

focused  on  the  microstructural  developments  in  metals  during  repeated 
5-9 

stressing.  During  the  1960s,  many  studies  were  performed  to  relate  the 

fatigue-induced  defect  structure  with  the  slip  morphology  exhibited  on 

...  ,  10-14 

the  surface. 


*A  complete  listing  of  references  is  given  on  page  29. 


1 


dtfftiiHiyiiftlifllfl 


While  surface  effects  were  recognized  as  a  controlling  factor  in 
fatigue  performance,  early  X-ray  methods  were  incapable  of  predicting 
the  impending  onset  of  fatigue  failure  or  even  the  span  of  the  fatigue 
life.^’^  Line  broadening  was  observed  after  cycling  over  a  small  frac¬ 
tion  of  the  total  fatigue  life,  but  then  it  remained  virtually  unaltered 

both  in  extent  and  intensity  throughout  the  remainder  of  the  life. 

17  18 

A  recent  application  of  a  special  X-ray  method  *  to  cyclically 
stressed  2024  A1  enabled  the  prediction  of  both  the  fatigue  life  and 
failure  of  the  aluminum  with  a  considerable  degree  of  accuracy.  This 
nondestructive  method  for  analysis  of  the  fatigue-induced  defect  structure 
is  based  on  the  principle  of  X-ray  double-crystal  dif fractometry  and  em¬ 
ploys  X-ray  topography  to  afford  a  visualization  of  the  defect  configura¬ 
tion.  The  polycrystalline  specimen  is  irradiated  with  a  crystal- 
monochromated  beam;  and  each  reflecting  grain  is  considered  to  function 
independently  as  the  test  crystal  of  a  double-crystal  diffractometer. 
Depending  on  the  perfection  of  the  grains,  the  specimen  is  rotated  in 
intervals  of  seconds  or  minutes  of  arc:  and  the  spot  reflections,  recorded 
along  the  Debye  arcs  of  a  cylindrical  film  for  each  discrete  specimen 
rotation,  are  separated  by  film  shifts.  This  multiple-exposure  technique 
gives  rise  to  an  array  of  spots  for  each  reflecting  grain.  These  arrays 
of  spots,  with  their  intensity  dependent  on  specimen  rotation,  represent 
X-ray  rocking  curves  of  the  reflecting  grains.  Thus,  if  the  grains  con¬ 
tain  a  substructure,  the  intensity  distribution  of  the  arrays  of  diffrac¬ 
tion  spots  will  be  multlpeaked  not  only  along  the  horizontal,  but  also 
along  the  azimuthal  elevation,  From  the  angle  subtending  successive  peaks 
of  the  rocking  curve,  the  excess  dislocation  density  between  subgrains 
can  be  determined;  while  the  excess  dislocation  density  within  the  sub- 
grain  can  be  obtained  from  the  spread  of  the  subpeak  curve.  From  the 

width  3,  at  half  the  maximum  intensity,  the  excess  dislocation  density  of 

17  18 

the  entire  grain  is  determined,  * 

By  analyzing  various  Chkg.)  reflections,  a  representative  statistical 
parameter  3  of  the  defect  structure  of  the  grain  population  is  obtained. 
Furthermore,  by  taking  Berg-Barrett  reflection  topographs  and  performing 


2 


a  spatial  tracing  at  the  reflections  to  the  spot  reflections  of  the 

rocking  curve,  the  analyzed  rocking  curve  can  be  correlated  to  the  grain 

18 

topography  on  the  specimen. 

The  double-crystal  diffractometer  arrangement  has  been  successfully 
employed  to  study  single  crystals,  pure  polycrystalline  metals,  and  multi¬ 
phase  alloys  with  ease.  The  method  offers  good  resolution  of  the  local 
effects  in  individual  grains,  but  also  provides  information  from  a  suffi¬ 
cient  number  of  grains  to  permit  accurate  statistical  analysis  of  an 
average  or  mean  microstructural  response.  The  film-recorded  reflections 
represent  an  imaging  of  the  subdomain  structure  which  can  also  be  inter¬ 
preted  quantitatively  in  terms  of  the  excess  dislocation  density,  which 
is  known  to  be  a  factor  in  the  accumulation  of  fatigue  damage  and  the 
ultimate  initiation  of  fracture.  The  purpose  of  the  present  study  is  to 
exploit  these  aspects  of  the  X-ray  method  to  isolate  a  suitable  micro- 
structural  indicator  of  progressive  fatigue  damage,  and  thereby  contribute 
to  the  capability  for  accurate  determination  of  the  fatigue  life  and 
failure  prediction.  In  accomplishing  this  objective,  the  double-crystal 
diffractometer  is  particularly  useful  for  the  microstructural  analysis 
of  both  the  surface  and  the  bulk.  Analysis  of  these  regions  can  be  car¬ 
ried  out  independently,  in  a  stepwise  manner,  by  incremental  removal  of 
surface  layers.  In  addition,  nondestructive  in-depth  analysis  is  made 
possible  by  altering  the  depth  of  penetration  of  the  incident  radiation. 
This  potential  may  be  very  important  in  technological  applications  when 
service  components  are  to  be  examined. 

A  method  which  can  be  used  in  conjunction  with  the  double-crystal 
technique  described  above  is  the  Warren- Averbach  (WA)  method  of  Fourier 
analysis  of  X-ray  line  broadening.  This  technique  is  unique  in  that  the 
sizes  of  coherently  reflecting  domains  and  the  internal  microstrains  can 
be  determined  simultaneously.  Alternate  theories  of  line  broadening 
assume  that  broadening  of  Bragg  peaks  is  due  predominately  to  either  size 
or  strain,  but  not  both.  The  WA  method  is  unique  in  that  multiple 
sources  of  line  broadening  may  be  operating  simultaneously. 


Two  peaks  are  selected  for  the  WA  analysis,  usually  the  200  and  400 
Bragg  reflections.  A  counter  is  used  to  record  the  intensity  distribu¬ 
tion,  and  the  intensity  peaks  at  0.005-degree  increments  are  fed  into  a 
PDP  11-34  computer.  The  Rachinger  graphical  separation  of  the  Ka^  and  Ka^ 
double  profile  is  employed  (see  Appendix  A) ;  and  the  Kot^  corrected  profile 
of  both  annealed  and  worked  samples  is  used  in  the  analysis.  The  Stokes 
method  is  used  for  the  correction  of  the  instrumental  broadening. 

A  FORTRAN  computer  program  was  written  (see  Appendix  B)  for  calculat¬ 
ing  the  Stokes-corrected  Fourier  coefficients  of  a  broadened  profile. 
According  to  the  WA  analysis,  the  corrected  Fourier  coefficients  are 
related  to  the  particle  size  and  the  root-mean-square  microstrains  as 
follows: 


*»  *1  <ho>  -  »"  4  -  <  £l  >  h0 

a 


(i) 


where 


*  the  corrected  Fourier  coefficients 


A^  =  the  corrected  Fourier  coefficients  of  particle  size 

L  *  the  distance  normal  to  the  reflecting  planes 

2  2  2  2 
hj  -  h  +  k4  +  V 

a  =  the  lattice  parameter 
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By  plotting  £n  A^(h^)  versus  h^  (usually  for  two  reflections)  the 
size  effect  and  microstrains  can  be  separated.  The  intercepts  of  these 

g 

plots  provide  the  particle  size  coefficients  A^  for  various  L  and  from 
the  slopes  the  microstrains  can  be  calculated. 


PREVIOUS  WORK 

Several  reports  covering  previous  work  under  this  program  have  been 
19-21 

published  recently.  These  elucidate  the  preferential  work  hardening 

occurring  at  the  surface  of  fatigued  AA2024-T4  samples  and  the  associated 
excess  dislocation  densities  for  grains  near  the  surface.  It  was  found 
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that  an  initial  rapid  increase  in  3  comprised  the  first  20  to  25  percent 
of  the  life  followed  by  a  long  intermediate  period  featuring  very  little 
structural  change  and  a  final  rapid  enhancement  during  the  last  5  to  10 
percent  of  the  life,  coinciding  with  crack  initiation  and  growth.  It 
was  found  that  a  critical  excess  dislocation  density  exists  which  is 
associated  with  fatigue  fracture.  This  critical  value  was  independent 
of  the  stress  amplitude  and  varied  according  to  a  Petch- type  relationship 
for  different  grain  sizes  of  the  same  alloy.  An  in-depth  analysis,  car¬ 
ried  out  by  stepwise  removal  of  the  surface  layers  of  cycled  specimens, 
disclosed  a  plastic  response  for  the  grains  located  in  the  bulk  after 
about  5  percent  of  the  fatigue  life.  The  defect  structure  in  the  speci¬ 
men  core  developed  gradually  during  the  cycling.  The  excess  dislocation 
density  for  the  bulk  increased  almost  linearly  as  a  function  of  the 
fatigue  life,  with  a  terminal  value  at  failure  identical  to  the  critical 
excess  dislocation  density  for  surface  grains.  The  fatigue  process  was 
interpreted  as  a  rapid  work  hardening  of  the  surface  to  form  a  barrier  to 
dislocation  egression  and  rearrangement.  The  dynamic  interplay  between 
the  surface  barrier  and  the  eventual  plastic  response  activated  in  the 
bulk  leads  to  a  critical  defect  accumulation  at  the  surface  and  incipient 
cracking. 

When  the  fatigue  process  was  interrupted  prior  to  failure  and  the 
surface  layer  was  removed,  a  striking  recovery  phenomenon  was  observed 
throughout  the  specimen  cross  section  during  subsequent  cycling.  The 
bulk  defect  structure  was  thus  shown  to  be  extremely  unstable  in  the 
absence  of  the  restraining  influence  imposed  by  the  work-hardened  surface 
layer.  The  extension  of  the  fatigue  life  of  metals  by  judicious  surface 
removal  was  ascribed  primarily  to  the  elimination  of  the  surface  barrier, 
rather  than  to  the  removal  of  microcracks. 

The  fatigue  response  at  various  depths  from  the  surface  was  also 
investigated  nondestructively  by  employing  X-ray  radiation  with  differing 
penetration  capabilities.  The  excess  dislocation  density  of  grains 
located  up  to  300  micrometers  in  depth  was  examined  using  molybdenum 
radiation  and  was  found  to  vary  linearly  with  the  fatigue  life.  This 
steep,  linear  dependence,  in  conjunction  with  the  early  life  saturation 


behavior  of  surface-grain  densities  measured  with  copper  radiation,  pro¬ 
vided  a  new  criterion  for  predicting  accurately  the  fatigue  life  and 
failure. 

PRESENT  WORK 

The  striking  feature  of  the  previous  work  is  that  a  linear  relation¬ 
ship  exists  between  the  statistical  average  $  and  the  remaining  fatigue 
life.  This  enabled  accurate  prediction  of  the  useful  remaining  life  by 
a  simple  nondestructive  determination  of  B  using  Mo  radiation.  However, 
these  results  were  strictly  applicable  only  to  samples  tested  at  a  con¬ 
stant  stress  amplitude.  Although  tests  indicated  that  (3*  was  stress 
amplitude  independent,  it  was  not  clear  how  a  single  sample  would  behave 
under  more  complex  loading  histories.  Accordingly,  the  present  tests 
were  undertaken  to  determine  the  applicability  of  the  X-ray  double-crystal 
dif fractometry  technique  for  predicting  fatigue  failure  in  samples  with 
a  random  stress  loading  sequence. 

The  spectral  loading  sequence  chosen  is  shown  in  Figure  1,  Four 
stress  levels  were  selected  between  25  ksi  and  41  ksl,  corresponding 
approximately  to  the  fatigue  limit  and  the  proportional  limit,  respec¬ 
tively.  The  number  of  cycles  at  each  stress  level  was  determined  by 
setting  each  loading  step  to  20  percent  of  the  fatigue  life  at  each  stress 
level.  The  stress  versus  number  of  cycles  (S-N)  curve  for  such  a  deter¬ 
mination  is  shown  in  Figure  2  and  was  obtained  for  a  constant  stress 
amplitude  (R  =  -1) . 

After  each  loading  step,  the  sample  was  analyzed  using  the  double¬ 
crystal  diffractometer  with  Mo  radiation  to  determine  the  rocking  curve 
halfwidth.  The  value  of  B  obtained  after  each  step  was  then  compared  to 
the  calibration  curve  shown  in  Figure  3.  This  curve  was  determined  in 
the  previous  studies  on  constant  amplitude  fatigue  cycling.  At  the  end 
of  the  fourth  loading  cycle  (41  ksi) ,  B  was  determined  and  located  on  the 
calibration  curve.  As  shown  in  Figure  4,  the  predicted  remaining  life 
was  39  percent  based  on  the  X-ray  method.  Figure  4  also  shows  the  pre¬ 
dicted  remaining  fatigue  life  of  20  percent  based  on  the  conventional 
technique  using  the  S-N  curve.  The  actual  remaining  fatigue  life  was 


found  to  be  42.5  percent.  Therefore,  the  error  in  estimating  the  remain¬ 
ing  life  on  the  basis  of  the  usual  techniques  was  over  100  percent; 
while  the  error  of  the  X-ray  technique  was  less  than  10  percent. 

Analyses  using  line  broadening  techniques  and  conventional  diffrac- 
tometrry  have  also  been  performed.  After  each  stage  of  cycling,  a  con¬ 
ventional  diffractometer  trace  is  run  to  determine  the  peak  halfwidths 
for  the  111,  200,  222,  and  400  peaks;  and  corrections  are  made  to  take 

into  account  variations  in  the  initial  halfwidths.  If  a  Cauchy  line  pro- 

22 

file  is  assumed,  it  can  be  shown  that  the  line  broadening  is  related 
to  the  coherently  reflecting  domain  size  8,  and  the  strain  e  through  the 
equation 

dSQ  =  j  +  2eS  (2) 

where  S  *  2sin0/A 

dS  =  2cos0d8 /A 

0  ®  Bragg  angle  for  hkJt  reflection 
A  =  X-ray  wavelength 

The  quantity  2d0  Is  equivalent  to  B,  the  peak  halfwidth;  and,  by  a 
simple  substitution.  Equation  (2)  becomes 

Bcos0  =  -J-  +  4esin0  (3) 

By  a  linear  regression  analysis,  the  quantities  £  and  e  can  be  readily 
determined  from  the  experimentally  measured  $.  As  shown  in  Figure  5,  the 
relationship  between  Bcos0  and  sin0  is  indeed  linear  for  the  four  frac¬ 
tions  of  life  investigated  (N/N-,  =  0.22,  0.28,  0.42,  and  0.57).  The 

r 

slope  of  each  line  is  equal  to  4e  and  the  intercept  is  A/8,. 

Once  the  values  of  e  and  8, are  determined,  the  individual  contribu¬ 
tions  to  the  dislocation  density  pg  and  can  be  determined  through  the 
equations 


(4a) 
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7,68  g2 
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PA 


(4b) 
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(pep4> 


1/2 


(4c) 


The  quantities  pg,  p^,  and  p  derived  from  Figure  5  are  summarized  in 
Table  1  along  with  the  quantities  e  and  Jt.  Note  that  the  coherently 
diffracting  domain  size  remains  nearly  constant  over  the  life  interval 
0.22  to  0.57,  while  the  strain  steadily  increases. 


TABLE  1  -  LINE  BROADENING  ANALYSIS  BY 
CONVENTIONAL  DIFFRACTOMETER 


Fraction 

of 

Life 

-(N/Np) 

Strain 

-e 

Domain 

Size 

o 

(A) 

-a 

Dislocation 

Density 

(cm~2) 

Pe 

Dislocation 

Density 

(cm”2) 

PJt 

Root  Mean 
Dislocation 

(cm~2) 

(pep//2 

0.22 

0.00157 

319 

2.33  x  1010 

9.83  x  1010 

4.79  x  10 10 

0.28 

0.00322 

462 

9.77  x  1010 

4.69  x  1010 

6.77  x  1010 

0.42 

0.00406 

344 

15.50  x  1010 

8.45  x  1010 

11.44  x  1010 

0.57 

0.00436 

305 

17.86  x  1010 

10.80  x  10 10 

13.89  x  10 10 
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|  Both  quantities  pg  and  p  are  linearly  related  to  the  fraction  of 

life  N/Np,  as  shown  in  Figures  6  and  7.  The  best  fit  is  obtained  by 

utilizing  p,  which  Includes  variations  in  both  p  and  p„ .  The  error  be- 

_  ex 

tween  the  experimental  value  of  p  and  the  value  obtained  from  the  linear 
regression  is  in  no  case  worse  than  10  percent.  Thus,  these  types  of 
measurements  should  be  capable  of  being  used  to  accurately  predict  the 
remaining  life.  Subsequent  experiments  are  being  initiated  to  verify 
this  assumption. 
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SPECTRAL  LOADING  SEQUENCE 


Spectral  Loading  Sequence  for  AA2024-T4  Samples 


Fatigue  Curve  Determined  from  Constant  Stress  Amplitude  Tests 


HiaiMJIVH  30VM3AV  Q2103HH00 


O  O  CM  00  ^ 

<N  «" 


(£oi*)  e  $09  d 


Results  of  Conventional  Diffractometer  Study  for 
'arious  Fractions  of  Fatigue  Life,  N/N^ 


FRACTION  OF  FATIGUE  LIFE 

Correlation  Between  Strain  Dislocation  Density  p 
and  Fraction  of  Fatigue  Life 


AOZ 


i 
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FRACTION  OF  FATIGUE  LIFE 

-  Correlation  Between  Mean  Dislocation  Density 
and  Fraction  of  Fatigue  Life 


THERE  EXISTS  DIFFERENCES  BETWEEN  THE 

UPDATED  HARD  COPY  AND  THIS  VERSION  OF  THE  PROGRAH 

PROGRAN  RACHC 
IMPLICIT  INTEGER  (A-Z) 

INTEGER  TLU.NP, IDEL2 

REAL  CPS ( 500 ) . DEL . IDEL . Y I NC  , SUM . A , B , C » FRA 

REAL*8  FSPEC 
REAL«4  DEV 

L0GICAL*1  DSN< 10) »TEM(3> 

DATA  DEV/3RDX1/ 

DO  401=1.10 
DSN(I>='  ' 

GO  TO  16 
WRITE(5,20> 

DO  431=1.10 
DSN ( I )  =  '  ' 

FORMAT(/'  FILENAME  MUST  BE  6  CHARACTERS  WITH  3  CHARACTER 
*  EXTENSION.',/.'  EX!  •FILNAM.DAT’  (OMIT  QUOTES)'.///.'  RE-') 

WRITE  (5,1) 

FORMAT ('  ENTER  THE  FILE  NAME*.',*) 

RE AD (5, 2) <DSN( I), 1=1, 10) 

FORMAT (10A1) 

1  =  1 

IF  (DSN(I).EQ.'.')  GO  TO  12 
1  =  1  +  1 

IF  (I.EO.IO)  GO  TO  11 
GO  TO  3 

IF  (I.GT.7)  GO  TO  15 
TEM< 1 )=DSN( I+l ) 

TEM  <  2) =DSN< 1+2 ) 

TEM<  3 ) =DSN ( 1+3 ) 

DSN<  7  >  =TEM( 1 ) 

DSN(8)=TEM<2.» 

DSN<  9 ) =TEM<  3 ) 

IF  (I.E0.7)  GO  TO  11 
DO  99K=I,6 
DSN ( K ) = '  ' 

DSN< 10)='  ' 

J  =  I RAD50  <  9 » DSN , FSPEC ) 

I  =  IASIGNC9. DEV, FSPEC) 

TLU=5 

FORMAT < / , '  ENTER  THE  2  THETA  SEPERATION  BETWEEN  THE 
•ALPHA  1  AND 

FORMAT < / , '  ENTER  THE  RATIO  BETWEEN  THE  ALPHA  1  AND  ALPHA  2(',*> 
FORMAT (/»  '  INTEGRATED  AREA  CORRECTED  FOR  ALPHA  l.  ALPHA  2  =',F9.2) 
READ (9,97) NP 
FORMAT (13) 

D030I=1 »NP 
READ ( 9 » 96 ) CPS ( I ) 

F0RMAT<F12.6> 

CONTINUE 

A=0. 

B=0. 

D042I=1 » 10 
A=A+CPS(I>/10. 

B=B+CPS(NP+1-I)/10. 

P038I-1.NP 

CPS(I)=CPS(I)-(A+(B-A)*(I-l)/(NP-l. ) ) 

WRITE(TLU« 17) 

READ (TLU.6) DEL 
FORMAT (F12.6) 

WRITE (5,7) 

FORMAT (/,  '  ENTER  YINC1',*) 

READ<5,8) YINC 
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FORMAT (F12.6) 

IDEL-DEL/YINC 
DEL-DEL/YINC 
C-DEL-IDEL 
SUM-0.0 

IDEL2-INT  < I DEL ) 

DOlOOI-l. IDEL2 
SUM-SUM+CPS ( I > 

WRITE<TLU»37 > 

READ<TLU »9>FRA 
FORMAT  <F12. A) 

DOllOI-IDEL+1  rNF' 

CPS< I )=CPS< I ) - <CPS ( I-IDEL+1 >-<CPS< I-IDEL+1 )-CPS( I-IDEL ) ) *C )*FRA 
SUM-SUM+CPS! I ) 

WRITE!5»9R) 

FORMAT ( / . '  THE  PRINTED  DATA  SET  WILL  BE  ALLOWED  TO  BE  STORED  ON' 
*/»'  DISC  AFTER  THE  RESULTS  HAVE  BEEN  VERIFIED.') 

D0350I-1.NP 
WRITE! A. 95) CPS! I ) 

CONTINUE 
FORMAT ( F 1 2 • 6 ) 

WRITE ( TLU » 41 ) SUM 
END 


APPENDIX  B 

FORTRAN  PROGRAM  FOR  CALCULATING  THE  STOKES- 
CORRECTED  FOURIER  COEFFICIENTS 
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REAL*-*  DEV 
REAL*4  DEV2 
REAL*4  FILE<  2) 

REAL*4  FILE2(2) 

DATA  DEV/3RDX1/ 

DATA  DEV2/3RDX1/ 

DATA  F ILE/6RU0RKED . 3RDAT / 

DATA  F ILE2/6RANEALD. 3RDAT/ 

I=IASIGN(9. DEV. FILE) 

I=IASIGN(8f DEV2.FILE2) 

DEFINE  FILE  9(500.4.U.NREC> 

DEFINE  FILE  8(500.4. U.NREC2) 

REAL  XINTA ( 500 ) » FRA < 60 ) «FIA(60> 

REAL  XIN7U(500>»FRU<60).FIU(60> 

REAL  FCCR(60> .FCCK60) .LNFCR(60> 

REAL  SAVE (60) 

REAL  LACA.LACU.KV.MA 
INTEGER*4  FLAGS (60) 

IN  IEGER*4  YES. NO. IANS 

INTEGER  TLU.IPRAM<5> , TOF . R1W< 128 ) »R1A(I28>  »TLEU(32> . 

♦  TLEA <  32 ) . DSU . DSA . PPDATA ( 3 ) . HKL . RADA . RADU . DATE  <  4 ) . 

♦  DATE A <  4 ) , DS . RAGE . HKLU <  4 ) . HL ( 4 ) 

REAL*8  DATEU.TIME 

EOUI VALENCE  < R1 A < 1 ) , DSA ) . ( R1 A< 2 ) . TLEA > . < R1 A< 34 ) . HL ) » 

♦  <R1A(40) . RADA ) . ( R1 A ( 48 ) . TA  > . < R1 A< 50 > .LACA ) . ( R1 A ( 52 > . 

♦YAINC) . <R1A(54) .TTHA1 ) » ( R1 A  ( 56  >  .  NF'A ) . (R1A(62> . DATEA > . ( R1 A ( 66 > , 
♦WAVE ) 

EQUIVALENCE  <R1W<1 ) .DSU) , <R1W<2> . fLEU) . (R1U(34) . 

♦HKLU) . <R1U< 40) .RADU) . (R1U(42).NA> . ( R1U< 44 > . KV ) . 

♦  ( R1U<  48 )  .  TU )  .  < RIU<  50 ) . LACU ) .  <  R1U <  52 )  »  YINCU )  .  <  R1U< 54 ) .  TTHU1 ) 

♦. (R1U(56) .NPU) . (R1U(62) . DATEU) . (R1U<66> .UUAVE) 

DATA  YES/3HYES/ 

DATA  N0/3HN0  / 

DATA  PPDATA/2HPP . 2HDA . 2HTA/ 

KV=0 . 0 

DO  66681=1.60 
6668  FLAGS ( I )  = '  ' 

1  FORMAT ('  THE  FOLLOUING  REQUESTS  ARE  FOR  THE  ANNEALED  SAMPLE.') 

2  FORMAT < '  UHAT  IS  THE  UAVELENGTH  USED') 

5  FORMAT < '  ENTER  THE  MAXIMUM  AND  MINIMUM  TUO  THETA  VALUES' 

♦//'  MAX=_'.*) 

6  FORMAT <  '  MIN=„.  .*> 

9  FORMAT < '  R.C.F.C.=REAL  CORRECTED  FOURIER  COEF.'./. 

♦  '  I ,C.F.C.=IMAGINARY  CORRECTED  FOURIER  COEF.'./. 

♦  '  L.F.C.C, =LOG  OF  THE  REAL  CORRECTED  FOURIER  COEF.'./. 

♦  '  DI3P.  “DISPLACEMENT  BETUEEN  SAMPLED  CELLS'./. 

♦  '  R.U.F.C.=REAL  UNCORRECTED  FOURIER  COEF.'./. 

♦  '  I. U.F.C. “IMAGINARY  UNCORRECTED  FOURIER  COEF.'./. 

♦  '  R.F.C.A. “REAL  FOURIER  COEF.  OF  ANNEALED  SAMPLE'./. 

♦  '  I. F.C. A. “IMAGINARY  FOURIER  COEF.  OF  ANNEALED  SAMPLE'./. 

♦  '  L.U.F «C.R.=L(iG  OF  UN -NORMALIZED  REAL  FOURIER  COEF.'.///) 

10  FORMAT  < '  HARMONIC  R.C.F.C.  I.C.F.C.  L.F.C.C.  DISP. 

♦  R. U.F.C.  I. U.F.C.  R.F.C.A.  I. F.C. A.  L.U.F.C.R.') 

16  FORMAT (2X.F5.2»2X.F8.3.2X»F8.3»6X.F6.3»6X»F7.2»2X»F8.3» 

♦5X  »F8. 3.5X.F8.3.3X .F8.3.3X.F8.3.3X.AI) 

18  FORMAT (IX) 

19  FORMAT < '  ENTER  THE  MINIMUM  TUO  THETA  VALUE  UHICH  UILL  EVER  BE 

♦  './.'  USED  FOR  THE  BROADEST  PEAK  USING  THIS  RADIATION') 

20  FORMAT ( '  ENTER  THE  BRAGG  ANGLE  FOR  THE  PEAK') 

21  FORMAT ( '  AAA“'.F12.7> 

31  FORMAT ( '  DO  YOU  UISH  TO  ENTER  DATA  FOR  ANOTHER  UORKED  PEAK 

♦  OF  THIS  REFLECTION'?') 

32  FORMAT (A3) 

33  F ORMAT ( '  THE  FOLLOUING  DATA  REQUESIS  ARE  FOR  THE  UORKED  SAMPLE') 

34  FORMAT ( '  DO  YOU  UISH  TO  CHANGE  THE  HARMONIC  FACTOR  FROM  5?') 
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38  FORMAT  < '  HARMONIC  FACTOR-  ') 

AO  FORMAT < '  ENTER  THE  SAMPLES  IDENTIFICATION') 

41  FORMAT <32A2> 

47  FORMAT (//, '  PAGE  ' . I3./.32A2, /. '  DATA  WAS  TAKEN  '.A9»5X»A5> 

48  FORMAT (35X, '  INITIAL  '.11X.'  THICKNESS' ./. 15X, '  RADIATION' »UX» 

♦  '2  THETA'. 13X,'IN  CM. ' ./. 19X.A2. 15X,F6.2,14X»F8.5,//,15X» 

*'  MILLIAMPS'.UX. 'INCREMENT' .11X. 'ABS.  COEF . ' ./. 19X.F4. 1 . 17X, 
♦F4.2.12X.F7.2.//.15X. '  KV  '. 12X. ' POINTS './. 18X.F4. 1 . 15X 
#,13»//.15X. '  MONOCHROMETER' .7X» 'AAA' »17X. 'HARMONIC' »/» 15X, 

*'  BRAGG  ANGLE', 30X. 'FACTOR' ,/, 18X,F6.2»8X,F12.6, 14X»F6.3»//> 

49  FORMAT <'  IS  THERE  ANOTHER  ANNEALED  PEAK  YOU  WISH  TO  USE 

•  FOR  A  STANDARD  ?') 

59  LACA=0. 

TA=0. 

PAGE=1 
WRIT  E  <  5. 3 ) 

3  FORMAT ( '  ENTER  THE  DATES',*) 

READ<5»4)DATEU 

4  FORMAT <A9> 

WRITERS, 7) 

7  FORMAT < '  ENTER  THE  TIMES',*) 

READ( 5,8) TIME 

8  FORMAT <A5) 

WRITE<5. 19) 

READ<5« 1111 >TTLQW 

1111  FORMAT (F12.6 ) 

TL0W«TTL0W*3. 1416/360. 

WRITE (5,20) 

READ  <5,1111) BRAGG 
BRAGG=BRAGG*3. 1416/180 
WRITE* 5, 1 ) 

86  WRITE(S,40) 

READ  (5, 41 )  <  TLEA<  I), 1  =  1, 32) 

WRITEC5.2) 

RE AD  <  5 , 1 1 1 1 ) WAVE 

PAUSE  'MOUNT  THE  DISC  CONTAINING  THE  ANNEALED  SAMPLE' 

READ* 8' 1 >NPA 
D023I=2.NPA+1 
READ<8' I)XINTA(I-l) 

23  CONTINUE 

24  WRITE (5, 5) 

READ(5,11 11 >TTHA2 
WRITE <5, 6) 

READ  <5,1111)TTHA1 
YAINC=<TTHA2-TTHA1 >/<NPA-l > 

350  XMXA-0.00 

001001=1, NPA-2 

DINT=XINTA< I )+XINTA< 1+1 )+XINTA< 1+2) 

IF ( XMXA.GT .DINT )GOT0100 
XMXA-DINT 
IMXA-.I  +  l 
100  CONTINUE 

WRI TE< 5, 3333) WAVE, BRAGG. TLOW 

3333  FORMAT < '  WAVE-' ,F12.6, '  BRAGG- ' ,F12.6» ' TLOW-' .F12.6 > 

AAA-WAVE/<4.0*<SIN(BRA0G)-SIN(Tl0W) ) > 

WRITE<5, 21 ) AAA 
HF-5. 

90  WRITE(5,34) 

READ(5,32)I ANS 
IF ( IANS.EQ.N0>G0r092 
IF  < IANS.NE • YES)G0T090 
WRITE(S,38) 

READ (5, 1111 )HF 
92  WRITE  <5.451) 

451  FORMAT < '  WHAT  RADIATION  WAS  USED, 

♦(ENTER  2-LETTER  CODF  FOR  TARGET  ELEMENT).',/,'  ?',*) 


23 


RE  AD  ( S  >  452  )  RADW 
452  FORMAT <A2> 

CALL  BKGND <  XINTA .NPA ) 

CALL  LRNTZ< WAVE* TTHA1.NPA. XINTA. YAINC.TLU.TA.LACA. ALPHA) 
CALL  FRC0FUMXA.XINTA.NPA.FRA.FIA.AAA.YA1NC.TTHA1.TTHA2. 
♦WAVE.TLU.HF) 

WRITE(5>33) 

370  LACW=0.0 

TW=0.0 

35  WRITE<5.40) 

READ<5.41) <TLEW(K) .K=l»32> 

WRITE<5. 2) 

READ <5. 111 1 > WWAVE 
00341=1 .500 
34  XINTW< I ) =0 . 0 

PAUSE  'MOUNT  THE  DISC  CONTAINING  THE  WORKED  SAMPLE ' 

READ<9'1)NPW 

0025 1*2  »  NPW+ 1 

READ<9' 1 ) XINTW ( 1-1 ) 

25  CONTINUE 

37  WRITE<5»5> 

READ<5>  111DTTHU2 
WRITE <5, 6) 

READ  <  5 . 1 1 1 1 ) TTHW1 
Y INCUS < T THU2-TTHU1 ) / (NPU-l ) 

102  XMXU-O. 0 

001051=1 .NPW-2 

UINT=XINTU(I )+XINTW< 1+1 >+XINTW( I+2> 

IF (XMXU.GT . DINT ) GOTO 105 
XMXU=PINT 
IMXU=1+1 
105  CONTINUE 

CALL  PKGND<XINTW.NPW> 

CALL  LRNTZ < WWAVE . TTHU1 » NPU .XINTW. YINCU. TLU. TW* LACW . ALPHA ) 
CALL  FRCOF < IMXW . XINTW . NPU. FRW.FIW . AAA . Y INCU . TTHU1 . TTHU2 . 
•WWAVE.TLU.HF) 

WRITE! 4. 47)PAGE,TLEW.DATEW. TIME 

WRI TE  <  4 . 48 ) RADW . TTLOW . TW . MA  »  Y INCW . LACW . KV . NPU . ALPHA , AAA . HF 
BAG=BRAGG* 1 80 . /3 . 1 4 1 4 
40  WRITE < 4.9) 

WRITE<4. 10) 

SLPMX=0. 

001701=1.40 

IF  <f'RA<  I )  .EQ.O,  .AND.FIAU)  .EQ.O.  >G0T0150 
FCCR(I)=<FRW<I)*FRA<I)+FIW<I)*FIA(I))/<FRA<I)**2+FIA(I) 

•«*2> 

FCCI<I)=<FRU<I)#FIA<I>+FIWU>*FRA(I) >/ <FRA< I >**2+FIA< I > 
***2> 

GOTO 155 

150  FCCR< I >=0. 

FCCI < I )=0. 

155  IF<I.LT.7>N=2#<I-1> 

IF<I.GT.6.AND.I.LT.12)N*5*<I-4> 

IF<I.GT.ll)N=lO*<I-8> 

ZZ=N/<AAA*HF> 

OIS«AAA*ZZ 

IF< I .EG. 1 )G0T0170 

DDID-DIS-DISO 

SLP=<FCCR< 1-1 )-FCCR< I ) )/DDIS 

IF<SLP.LT.SLPMX)G0T0170 

SLPMX-SLP 

AO=FCCR< I ) +DIS*SLPMX 
170  DISO-DIS 

0044491=1.40 

IF  <FCCR<I).LE.O. )  SAVE ( I ) =FCCR  < I ) 

IF  <FCCR(I).LE.O. )  FLAGS! !)■'*' 


IF  <FCCR< I ) .LE.O. )  GO  TO  6669 
SAVE ( I ) -ALOG < FCCR (I >  > 

6669  CONTINUE 

001751=1.60 
FCCR< I >=FCCR< I >/A0 
FCCI < I )=FCCI <  I ) /AO 
LNFCR<I>=-11. 51293 

IF ( FCCR  < I > . GT . 0 . 00001 ) LNFCR < I > =ALOG<  FCCR  < I > ) 

172  IF!I.LT.7)N=2*<I-1> 

IF!I.GT.6,AND.I.LT.12>N=5*< 1-4) 

IF!I.GT.ll)N=10*<I-8> 

Z2=N/ <  AAA*HF ) 

DIS=AAA*ZZ 

WRITE (6. 16)ZZ,FCCR< I ) .FCCI (I) . LNFCR < I > .DIS.FRW! I > . 

♦FIW< I ) »FRA< I ) »FIA< I ) .SAVE! I ) .FLAGS! I > 

IF!I/5*5.EQ.I)WRITEC6,18> 

IF! <1-25 >/45*45.NE. 1-25) G0T0175 
PAGE=PAGE+1 

WRI TE < 6 . 47 ) PAGE . TLEW. DATEW. TIME 
WRITE (6.10) 

175  CONTINUE 

210  URITE<5.31> 

READ <5. 32) IANS 
t  F  < I ANS . EQ . YES ) G0T0370 
I F  < I ANS . NE . NO  >  GO T021 0 
215  WRITE  <  5 » 49 ) 

RE AD <5. 32) IANS 
IF ( I ANS. EQ. YES) G0T059 
IF < I ANS.NE ■ NO) GQTQ215 
220  STOP 

END 

SllBROUTINEFRCOF< IMAX.XINT.NP.FR, FI. AAA, YINC.TTH1. 
♦TTH2.WAVE.TLU.HF) 

WRITE (5. 4444 ) 

4444  FORMAT!/,'  *FRCOF'> 

REAL  XSYM<251 ) ,XASYM<251 > ,XINT<500> ,FR<60>  »FI !60> 

1  FORMAT!'  YNORM=' .F10.6) 

5  FORMAT! '  THE  NUMBER  OF  POINTS  IN  RECIPROCAL  SPACE', 15, 

#/,'  HAS  EXCEEDED  THE  PERMITTED  VALUE  OF  251.  AN  ERROR  WILL 
RESULT.  CHANGE  THE  DIMENSIONS  OF  XSYM  AND  XASYM ' > 

8  FORMAT!'  THE  MAXIMUM  VALUE  OF  THE  CALCULATED  RECIPROCAL 
♦SPACE  COORDINATE' ,F8. 3,/, '  HAS  EXCEEDED  0.5.  THE  FOURIER 
♦  COEF.S  WHICH  FOLLOW  CONTAIN  THIS  ERROR.') 

RINC=YINC*3. 1416/<2. 0*180. ) 

TMX=  <  <TTH1+YINC*< IMAX-1 ) >/2. >*3.1416/180. 

DELS=2 . *AAA* < S IN <  TMX+R INC ) -S IN < TMX >  > /WAVE 

SI =2 . *AAA*S IN  <  TMX ) /WAVE 

S2=Si 

DO  91=1,251 
XASYM! I >=0.0 

9  XSYM!  I >=0.0 

XSYM! 1 ) =2, *X I NT! I MAX) 

XASYM! 1 )=0 .0 

MM=2.*AAA* (SIN!TMX)-SIN!TTHl*3.1416/<2 . *180. > ) > / < DELS*WA VE > + 1 
MMM=2.*AAA* (SIN! TTH2*3 .1416/ <2. *180. ) >-SIN! TMX ) ) /<  DELS*WAVE ) +1 
IFIMMM.GT .MM)MM=MMM 
N=1 

10  IF!MM.LE.251 )G0T015 
WRITE <5, 5) MM 

15  D0110I=1»MM-1 

FLG1-0. 

FLG2=0. 

S1*S1-DELS 
S2*S2+DELS 
T3*Sl *WAVE / <  2 . *AA A ) 

T5=S0RT!ABS< 1-T3*T3> ) 
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T1»ATAN!T3/T5> 

T4= ! S24WA0E/ ( 2 . * AAA ) ) 

TA=SQRT!1-T4*T4> 

T2=ATAN!T4/TA> 

M-NP-IMAX 

IF<IMAX.GT.M)M«IMAX-1 

XHIGH-0.0 

XLOU-O.O 

D0100J-N.H+3 

TH1=TMX- J*RINC 

TH2=TMX+J*RINC 

IF<TH1.GT.  T  i.Qft.FLGl.HE.  0.  >GOT050 

JN=J 

FLG1-1.0 

IF!IMAX-J.GT.0)G0T020 

XX1=0.0 

fioroaa 

20  XX1=XINT  < IMAX-J) 

IF ( XXI « LT . 0. )XX1=0. 

25  IF  < IMAX-J+1 .GT  <0)G0T030 

XLOU-O.O 
G0T050 

30  XX2=XINT ( IMAX-J+1 > 

IFIXX2.LT .0. )XX2=0. 
XL0U=XX1+!XX2-XX1)*!T1~TH1)/RINC 
50  IF( TH2.LT . T2.0R.FLG2.NE.0.OG0T090 

Jh=J 

FLG2-1.0 

IF< IMAX+J.LT  «NP+1 ) GOT 060 
YY 1=0.0 
GO T 065 

AO  YY1=XINT  < IhAX+J) 

IF!YY1.LT.O.)YY1=0. 

AS  IF  <  IMAX+J-1  .LT  .NF-+1  )GOT070 

XHIGH=0 . 0 
G0T090 

70  YY2=XINT  < IMAX+J-1 ) 

IF ( YY2.LT , 0. )YY2=0. 

XHIGH=YYl  +  <  YY2-YY1 >*! TH2-T21/RINC 
90  IF (FLU1 . NE . 1 . .OR . FLG2.NE . 1 . ) GOTO 100 

G0T0102 

100  CONTINUE 

102  XSYM! 1  +  1 ) =XHIGH+XLOU 

XASYM! 1+1 )=XHIGH-XLOW 
P=I 

RAT=XSYM!I  +  l)/C0S!3.14159265*DELS*P*HAR/!0.5>  > 
N»JN 

IFIJM.LT. JN)N=JM 
IF ( I .EQ. 1 )G0TQ110 
KI=I-1 

110  CONTINUE 

120  001301=1 »A0 

FR( I >=0.0 
130  FI(I)=0.0 

H=!KI-1)*DELS 
IF!Q.GT.0.5)WRITE!5.8)G 
WRITE !5» 5555) 

5555  FORMAT!//. '  A  LENTHY  CALCULATION  HAS  STARTED. 
♦'./»'  PROGRESS  IS  DOCUMENTED'.//.'  NOTES 
♦EXECUTION  TERMINATES  AT  O') 

ITEST-60 
D0150I-1 > AO 
WRITE<5.555A) I TEST 

5556  FORMAT  < '  -.12.*) 

IF < I  TEST .EQ. 30) WRITE <5. 5557) 

5557  FORMAT!'  ') 


ITEST-ITEST-1 

IF<I.LT.7)N-2*<I-1> 

IF<I.GT.A.AND.I.LT.12)N-5*<I-4> 

IF<I.GT.U)N-10*<I-8> 

Z-N/<AAA*HF) 

D0145J-1*KI 
S-DELS*< J-l) 

B-.AAAAAAA7 

IF(J/2*2.EQ.J)B-1. 3333333 
IF  <  J  .EQ .  1  >  B= . 33333333 

FR<  I ) -FR< I ) +XSYM  <  J ) *COS<  3 • 141 A*Z*S/0. 5 )*DELS*B 
145  FI<I)=FI<I)+XASYM<J)*SIN<3.1416*Z*S/0.5)*DELS*B 

IF ( I . EQ • 1 ) YNORM-FR  < 1 ) 

FR  <  I ) =FR  <  I ) /YNORM 
150  FI<I)«FI<I)/YNORM 

URITE<5» 1 > YNORM 
RETURN 
END 

SUBROUT I NEBKGND (XINT  »NP) 

URITE<5, 1 > 

I  FORMAT  < '  *BKGND ' ) 

REAL  XINT <500) 

BKGMN-0.0 

BKGMX-0.0 

D0120I=1,9 

BKGMN=BKGMN+X I NT  < I ) /9 . 0 
120  BKGMX-BKGMX+XINT<NP+l-I>/9.0 

BGSLF  =  <  BKGMX-BKGMN ) / <  NP- 10 > 

D0130I-1.NP 

XINT  < I ) =XINT  < I ) - < BKGMN+ BGSLP* < I - 5 ) > 

130  IF<I.LT.5.0R.I.GT.NP-5)XINT<I>-0.0 

RETURN 
END 

SUBROUTINELRNTZ< WAVE, TTHl.NP, XINT, YINC.TLU.T, LAC, ALPHA) 
URITE<5,11> 

II  FORMAT </,'  *LRNTZ' ) 

REAL  LAC, XINT <500) 

INTEGER44  YES, NO, IANS 
DATA  YES/3HYES/ 

DATA  N0/3HN0  / 

1  FORMAT < '  IS  THE  DATA  FROM  A  THIN  FILM?') 

2  FORMAT < '  ENTER  THE  LINEAR  ABSORPTION  CQEF.',/,'  L.A.C.- 

3  FORMAT < '  WHAT  IS  THE  FILM  THICKNESS  IN  CENTIMETERS' 

*./.'  T=_ ' ) 

4  FORMAT ('  NAS  A  MONOCHROMETER  USED?') 

5  FORMAT < '  ENTER  THE  BRAGG  ANGLE  OF  THE  MONOCHROMETER.',/ 

•'  ALPHA-,',*) 

6  FORMAT < A3) 

IF<  T . NE • 0)GOT030 
10  MRITE<5, 1 > 

READ <5, 6) IANS 
IF  < IANS.EQ.N0)G0T020 
IF< IANS . NE . YES ) GOTO 10 
WRITE<5,2> 

READ <5, 11 1 1 >LAC 
1111  FORMAT  <F12«A> 

WRITE<5,3> 

READ<5, 1111 )T 
G0T030 

20  LAC-4 

T-4 

30  ALPHA-0.0 

IANS- YES 

35  URITE<5,4> 

READ<5, A) IANS 

IF ( IAN8.E0.N0)G0T050 
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IF < IANS . NE . YES ) GOT030 
MRITE<5.5) 

RfcAIK5.ll  11 ) ALPHA 

SO  CSASQ- ( COS <  2* ALPHA* 3 . 1416/180 . 0  >  > **2 

D040I-1 »NP 
N-I-l 

TT«<TTH1+N*YINC>*3. 1416/180.0 
XZ*-2.0*LAC*T/SIN<TT/2.0> 

Tl“SIN<TT/2. )*SIN<  TT/2. > 

TS=COS(TT) 

T3“TS*TS 
1 2=C0S<  TT/2. ) 

T4=EXP<XZ> 

XINT<I)=XINT<I ) *SIN ( TT  > *<  SIN( TT/2 . 0 )/COS ( TT/2.0) >*( ( 

•  1 . 0+C0S<  TT )**2. 0 ) / ( 1 . 0+CSASQ4C0S ( TT ) **2) )/< 1 .0-EXP (-2.0 
**LAC*T/SIN<TT/2) > ) 

40  CONTINUE 

RETURN 
END 
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